Site-selective Mott transition in rare earth nickelates 
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A combination of density functional and dynamical mean held theory calculations are used to 
show that the remarkable metal-insulator transition in the rare earth nickelate perovskites arise 
from a site-selective Mott phase, in which the d-electrons on a half of the Ni ions are localized 
to form a fluctuating moment while the d-electrons on other Ni ions form a singlet with holes on 
the surrounding oxygen ions. The calculation reproduces key features observed in the nickelate 
materials, including an insulating gap in the paramagnetic state, a strong variation of static mag- 
netic moments among Ni sites and an absence of "charge order". A connection between structure 
and insulating behavior is documented. The site-selective Mott transition may be a more broadly 
applicable concept in the description of correlated materials. 
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Rare earth nickelates ReNi03 with Re=Sm, Eu, Y or 
Lu undergo a remarkable paramagnetic metal to para- 
magnetic insulator transition as the temperature is de- 
creased below a temperature Tmi ~ 400-600K The 
insulating state is characterized by a two-sublattice sym- 
metry breaking, with the Ni on one sublattice having a 
decreased Ni-0 bond length and the Ni on the other hav- 
ing an increased mean Ni-0 bond length. While the ma- 
terials arc sometimes classified as Mott/charge-transfcr 
materials [2 4], the Mott picture docs not account for 
the association of insulating behavior with bond dispro- 
portionation. Charge ordering has been extensively dis- 
cussed P, but charge ordering in the naive sense of 
a change in the Ni valence between the two sublattices is 
expected to be suppressed by the large repulsive d-d in- 
teraction U on the Ni site and appears to be ruled out by 
very recent soft X-ray resonant diffraction data An- 
tifcrromagnctic order can drive an insulating state and 
induce a disproportionation as a second-order effect [Io| ; 
while this mechanism may be relevant to NdNiC>3 , where 
the magnetic and metal-insulator transitions coincide, it 
fails to account for the paramagnetic insulating state ob- 
served in the Sm, Eu Y and Lu compounds. The origin 
of the paramagnetic insulating state has heretofore re- 
mained mysterious. 

Here we present density functional plus dynamical 
mean field theory (DFT+DMFT) calculations which 
show that the bond-length disproportionation and asso- 
ciated insulating behavior are signatures of a novel corre- 
lation effect in which the Ni d electrons on one sublattice 
(Nil, long bonds) become effectively decoupled from the 
surrounding lattice while the d electrons on the other 
sublattice (M2 , short bonds) bind with holes on the oxy- 
gen sites to form a singlet state. We term this state a 
site-selective Mott insulator . In the site-selective Mott 
state the singlet formation energy dictates the magni- 
tude of the insulating gap, while the valence difference 
between the two Ni sites, of course nonzero by symme- 
try, is very small and is not relevant to the physics. 



Our DFT+DMFT calculations are performed using 
the Vienna ab-initio simulation package [ll|, Ell to ob- 
tain bands from which localized Ni d orbitals (DMFT 
basis set) and O p orbitals are constructed using maxi- 
mally localized Wannier functions [l3| defined over the 
full ~40eV range spanned by the p-d band complex. In 
most of our calculations we used the LUMO3 crystal 
structure with the P2\/n space group and lattice param- 
eters taken from a fully relaxed DFT+U calculation of 
the magnetically ordered state; these are in close agree- 
ment with the experimentally measured parameters [l[ 
(See Table[T| for the antiferromagnetic and paramagnetic 
insulating states. In a few comparison calculations we 
used the Pbnm structure of the high-T phase of LUMO3. 
The LuNiOa crystal structure is such that e g and t 2g are 
no longer irreducible representations of the point group 
of either Ni site. However, for each Ni the oxygen cage 
is sufficiently close to cubic that one can still choose a 
coordinate system where there are a set of three and a 
set of two cubic harmonics which are closely grouped in 
energy, and the off-diagonal elements are relatively small. 
We determine the local coordinates for each d Wannier 
function such that the sum of the square of off-diagonal 
terms in the d manifolds on each site is minimized. We 
refer to the approximately twofold (threefold) degenerate 
states as e g (t 2g ) respectively. 

We treat the filled t 2g orbitals with a static Hartrec- 
Fock approximation (shown by recent work to be ade- 
quate except for fine details of the photoemission spec- 
trum [3]) while correlations in the Ni e g manifold are 
treated within single-site DMFT including the full ro- 
tationally invariant Slater-Kanamori interactions using 
t7=5eV and J=leV. The DMFT impurity problem is 
solved u sing th e continuous time quantum Monte-Carlo 
method [ljJH . An important issue in the DFT+DMFT 
procedure is the 'double counting correction' which ac- 
counts for the part of U already included in the under- 
lying DFT calculation and plays an important physical 
role by setting the mean energy difference between d and 
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Figure 1. (Color online) Momentum- integrated spectral 
function A(iS) computed for P2i/n LuNiC>3, projected onto 
the local Wannier basis corresponding to the Ni e 9 states for 
Nil (dashed, red) and N12 (solid, blue) and displayed in the 
near fermi-surface frequency regime, (a) DFT spectral func- 
tion, (b) DFT+DMFT spectral function obtained at (7=5eV, 
J=leV, 7V d =8.0, and T=116K. 



Figure 2. (Color online) Local magnetic susceptibility x(w = 
0) of LuNi0 3 calculated via DFT+DMFT as a function of 
temperature in the paramagnetic state using the low T P2\ jn 
structure with mean bond length difference of 0.104A. Red 
circles, Nil (large Ni-0 bond length). Blue square, N12 (small 
bond length). Inset: static magnetic moments n-f — 714. from 
DFT+DMFT calculations. 



p bands and thereby controlling the d occupancy. There 
is 110 exact procedure for determining the double count- 
ing correction; we therefore consider a range of values, 
parametrized by the resulting d-occupancy and chosen 
such that the d-valence within DFT+DMFT is not sub- 
stantially different from the valence found in the density 
functional calculations. We find that changes to U / J or 
the double counting correction shift the phase boundary 
slightly but do not change the main results. 

Fig. [T] compares the d-electron spectral function for 
P2i/n LUN1O3 obtained from a density functional calcu- 
lation (Fig. UK) to that obtained from our DFT+DMFT 
calculation (Fig. [lb). We see that while density func- 
tional theory predicts that the system is metallic (no gap 
at the chemical potential) despite the lattice distortion, 
the correlation effects captured by DMFT drive the sys- 
tem into an insulating state with a gap of ~ 200meV 
(the precise gap value depends on the double counting 
prescription and the value chosen for J). Differences in 
the near fermi-surface spectral functions for the two sites 
are evident in both cases, with the Nil site having larger 
spectral weight in the low energy regime below the fermi 
level and the N12 site having larger spectral weight above. 
This qualitative difference was interpreted by Ref. Q as 
evidence of charge ordering. It is important to note, how- 
ever, that the charge density is defined in terms of an in- 
tegral of the spectral function over the entire frequency 
regime. We find that the total d-charge (Wannier basis) 
is 8.22 ± 0.08 for the DFT spectra and 8.0 ± 0.06 for the 
DFT+DMFT spectra in Fig.Q] (here the + corresponds 



to Nil and — to N12); our DFT+U calculations (atomic 
orbital basis, see TableflJ give an even smaller difference. 
Thus we conclude that the difference in charge between 
the two sublattices, while of course nonzero by symme- 
try, is not physically important and in particular is not 
the cause of the insulating behavior. 

Fig. [2] shows the temperature dependence of the local 
magnetic susceptibility. A dramatic difference is visible 
between the two sites. The local susceptibility of the Nil 
site has a clear Curie (~ 1 /T) behavior indicating well de- 
fined and long-lived but thermally fluctuating magnetic 
moments, as expected in a paramagnetic Mott insulator 
(in the single-site DMFT approximation). By contrast, 
the local susceptibility of the Ni 2 site has a negligible 
temperature dependence, indicating a lack of a local mo- 
ment on this site. For this reason, we identify the origin 
of the insulating state as a site-selective Mott transition . 

Further insight into the physics of the insulating state 
is obtained from the self energy, shown in Fig. EI 

In the high temperature Pbnm structure the Nil and 
N12 atoms are equivalent. The imaginary part of S(w) 
(dashed lines, black on-line) shows the quadratic varia- 
tion expected of a Fermi liquid while the real part (not 
shown here) has a linear frequency dependence at low en- 
ergy indicating that the system is moderately correlated 
with a mass renormalization of 2.46, consistent with op- 
tical data on paramagnetic metallic LaNi03 18[. When 
the size of the disproportionation is increased to a mean 
Ni-0 bond length difference Sa = 0.104A, comparable 
to the distortion in the experimental structure, the so- 
lution becomes insulating and the self energy changes. 
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Figure 3. (Color online) The imaginary part of the self 
energies ImS(iii) of LuNi03 for the Nil site (panel (a), red) 
and Ni 2 site (panel (b), blue) calculated via DFT+DMFT 
in the paramagnetic state using the low T P2\/n structure 
with mean bond length difference of 0.1044 and analytically 
continued following Ref. [13]. The self energies in the high 
T Pbnm structure with the mean bond length difference as 
0.0A (black dashed line) are also given for comparison. 



Figure 4. (Color online) The metal-insulator phase diagram 
of LuNiOs in the plane of bond disproportionation 5a and 
double counting (here parametrized by mean d-ocupancy Nd), 
obtained by DFT+DMFT (red circle) and DFT+U (green 
diamonds). The lattice disproportionation and mean d occu- 
pancy obtained from our DFT+U calculation are shown as 
the vertical and horizontal dashed lines respectively. 



For the Nil site, the imaginary part of the self energy 
displays a pole at zero energy, as expected for a Mott in- 
sulator. However, the self energy of the N12 site displays 
a gap, with no pole, as expected in a singlet state with 
the quantum numbers of a trivial insulator. 

The different behavior of the two sites can be under- 
stood by considering the limit of extreme disproportion- 
ation, in which the Nil site is completely decoupled from 
the surrounding oxygen ions while the N12 site is strongly 
coupled. The mean Ni valence found in our and other 
DFT calculations is approximately d 8 (two electrons in 
the e g orbitals on each Ni and one hole per Ni in the oxy- 
gen orbitals). On each Ni site, the Hunds coupling results 
in a high spin d 8 (total spin 5 = 1) configuration. On 
the long-bond site a spin 1 local moment results; on the 
short-bond M2 site the Ni spin-1 is strongly coupled to 
the two holes on the O sites of the octahedron, forming a 
singlet, so the local susceptibility is T-independent, while 
the corresponding self energy has a gap (as in Kondo in- 
sulators). This model shows that an insulator can be ob- 
tained even though the variation in charge between the 
two Ni sites is negligible and the oxygen hole density is 
equally distributed over the oxygen network; hence this is 
not a charge ordering. A physical picture similar to what 
is discussed here was mentioned as one of several possibil- 
ities in a Hartree-Fock study by Mizokawa, Khomskii and 
Sawatzky [l9[, and has been discussed by Freeland [20j . 

To determine the relation between the lattice distor- 
tion and the insulating phase, we studied a sequence 
of crystal structures interpolating between the high T 



undistorted phase and the low-T distorted structure, 
parametrizing the structure by the mean bond dispro- 
portionation da. Fig. 2] shows a phase diagram in the 
space of 5a and the double counting correction (here 
parameterized by the d-occupancy N d ) computed using 
our DFT+DMFT procedure and, for comparison, via 
DFT+U which is (up to minor implementation differ- 
ences not relevant here) the Hartree-Fock approximation 
to DFT+DMFT. We used the VASP implementation of 
DFT+U, but with with the double counting correction 
modified from the standard VASP E%, = U(N d - 0.5) - 
J(N% — 0.5) (a is the spin index) to 

E° ic = U(N d - 0.5) - J(NS - 0.5) + E A (1) 

In the VASP DFT+U implementation, the d-occupancy 
is fixed by a charge self-consistency loop; variation of 
-Ea changes the final value, enabling computation of the 
phase diagram. The d-orbitals are defined slightly differ- 
ently in the two methods (via maximally localized Wan- 
nier functions in DFT+DMFT and via projection onto 
atomic d orbital functions defined within a sphere (radius 
set by the VASP default) around Ni atoms in DFT+U). 
As can be seen from Fig.f?]and by comparing the numbers 
in Table|I] to those quoted above, these differences are not 
relevant to the physics we discuss here: the basic size of 
the (i-occupancy, the difference in e?-occupancy between 
difference sites, and the variation with double counting 
correction are essentially the same. 

Fig. @] shows that in the case of DMFT, changing the 
double counting correction so as to increase the d- valence 
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Ni-0 bond[A] 


Mom [/is] 


N d 


M/I 


DFT 


1.94/1.94 
1.96/2.01 
2.01/1.96 


0.74/0.73 


8.21/8.20 


M 


DFT+U 


1.99/1.90 
2.02/1.94 
2.04/1.91 


1.52/0.57 


8.22/8.24 


I 


DFT+U 
(Jahn- Teller) 


1.92/1.92 
1.94/2.06 
2.06/1.94 


1.13/1.13 


8.22/8.22 


M 


Exp. 
(Pbnm) 


1.93/1.93 
1.96/2.00 
2.00/1.96 


N/A 


N/A 


M 


Exp. 

(P2i/n) 


1.98/1.89 
2.00/1.94 
2.03/1.92 


N/A 


N/A 


I 



Table I. The Ni-O bond length, magnetic moment, and d 
occupancy for LuNi03 obtained from fully relaxed DFT and 
DFT+U structural calculations. The Nd values are obtained 
by projecting the Kohn-Sham wavefunctions to atomic d or- 
bital functions defined within a sphere of radius 1.3A around 
Ni atoms. For the DFT+U calculations, U=5eV and J=leV 
are used. The Ni-0 bond lengths for both the high-T (Pbnm) 
and low-T (P2\/n) experimental structures [l| are also given 
for comparision. 



away from d s increases the distortion magnitude required 
to obtain the insulating phase. The DFT+U curve is 
much steeper, and favors the insulating state, as expected 
from the Hartree-Fock nature of this approximation. 

Finally, we consider the low temperature phase tran- 
sition from the paramagnetic insulator to the magnetic 
insulator (AFI) (both having the P2\/n structure). The 
inset of Fig. H shows the DFT+DMFT results for the 
temperature dependence of the magnetic moments on the 
two sublattices as the temperature is decreased below the 
Curie temperature. The low T limiting values are similar 
to experiment and to other calculations, but the magnetic 
transition temperature is grossly overestimated, presum- 
ably because of the neglect of spatial fluctuations and 

Additionally, the 



correlations in single-site DMFT [21 1 
ordered state is wrongly predicted to be ferromagnetic, 
as is also found in DFT+U calculations which will be 
analyzed further below. 



Although they cannot represent the paramagnetic in- 
sulating state, static mean field theories such as DFT, 
DFT+U and hybrid functional approaches may cap- 
ture some of the physics of the antiferromagnetic insu- 
lating ground state. Performing full structural relax- 
ations within spin polarized DFT shows metallic behavior 
with neither magnetism nor bond disproportionation (see 
Tabled]): the relaxed DFT result closely matches the high 
temperature Pbnm structure. This qualitative structural 
error is uncommon for DFT and signals the importance 



of correlations. We have also performed DFT+U calcu- 
lations; this approximation stabilizes the disproportion- 
ated structure, giving structural parameters very similar 
to experiment. The ground state of the DFT+U solution 
is found to be insulating and magnetic; the magnitude of 
the moments is substantially different between the two 
sites, with the Ni atom with an expanded octahedron 
having a moment of 1.52/is and the Ni atom with a con- 
tracted octahedron having a moment of 0.57/is, similar 
to the DFT+DMFT calculation and to what was found 
experimentally [B|,[22j and previously computed 0- How- 
ever, in both DFT+U and DFT+DMFT the magnetic 
ordering is found to be ferromagnetic; the experimen- 
tally observed antiferromagnetic state is metastable, be- 
ing some 30meV/Ni higher in energy in DFT+U. Fur- 
thermore, in the DFT+U metastable antiferromagnetic 
solution the smaller moment is strongly suppressed to ~ 
0.02/is, which is much smaller than the measured value. 
Very recently, variational self-interaction-corrected den- 
sity functional calculations [23[ in LaNiOa/LaAlOs su- 
perlattices were shown to produce the experimentally ob- 
served antiferromagnetic ground state within the nicke- 
late layer, though the smaller moment is again much too 
small relative to experiment. This suggests that non- 
local correlations may be necessary to properly describe 
the magnetic ordering, but we emphasize that the form 
of the ordered state is merely a low energy detail, clearly 
not responsible for the structural or metal-insulator tran- 
sitions. We also emphasize that in all of the methods and 
d orbital definitions we have explored, the d charge differ- 
ence between the Nii/Ni2 sites is found to be negligible: 
the insulating state should not be interpreted as charge- 
ordered. Finally, within DFT+U there is a metastable 
state with an equivalent Jahn- Teller distortion on each 
Ni site very close (~ 6meV/Ni) in energy It is not im- 
possible that this state could be realized under different 
circumstances. 

In conclusion, we have identified a new mechanism for 
a correlation-driven metal-insulator transition, in which 
a translation symmetry breaking lattice distortion leads 
to different spin physics (local moment vs singlet for- 
mation) on the inequivalent transition metal sites of the 
insulator. The mechanism is shown to account for the 
essential properties of the paramagnetic insulating state 
of LuNiC>3 and related materials. A key prediction of 
our work is a strong site dependence of the local mag- 
netic susceptibility, with one sublatticc exhibiting a Curie 
behavior and the other a temperature independent \. 
The transition we have identified may be related to the 
insulating behavior observed in few-layer LaNi03-based 
heterostructures 24, 25 1 where DFT+U calculations [25[ 
suggest that disproportionation can occur and drive a 
transition. 
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